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Abstract 

Mapping the magnetic exchange interactions from model Hamiltonian to density functional the- 
ory is a crucial step in multi-scale modeling calculations. Considering the usual magnetic force 
theorem but with arbitrary rotational angles of the spin moments, a spurious anisotropy in the 
standard mapping procedure is shown to occur provided by bilinear-like contributions of high order 
spin interactions. The evaluation of this anisotropy gives a hint on the strength of non-bilinear 
terms characterizing the system under investigation. 
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Multiscale modeling approaches are extremely important for describing huge magnetic 
systems, e.g. at the micrometer-scale which would be impossible with only density functional 
theory (DFT). In magnetism, usually the multiscale approach is performed after mapping 
the magnetic exchange interactions (MEI) of a classical Heisenberg model to the DFT coun- 
terparts. This is a crucial task which can lead to wrong results if not done carefully. The 
simple model is described by 



H = - 



9 J *i 



(1) 



where describes the pairwise (two-spin) MEI between spins at lattice sites i and j while 
ei(l,9,4>) defines the direction of the local moment Mj. Sometimes, higher order terms such 
as the four-spin or the biquadratic MEI are introduced in the previous Hamiltonian for a 
better mapping of the DFT results 

HQ. 

Once the MEI extracted, the investigation of magnetism of several type of systems can 

n n n n q 

be performed going from molecules|3], transition metals alloys|4j, [5] and surfaces[6|, |7|, di- 
luted magnetic semiconductors [§, [J, to clusters 10 13 1 and even for strongly correlated 



stytems 



141 ]. Thermodynamical properties are then easily accessible such as Curie tem- 



peratures, specific heat or magnetic excitation spectra and spin waves stiffness in multi- 
dimensional systems. 

An elegant method to extract the MEI is based on a Green function technique which 
has been derived 20 years ago by Lichtenstein and coworkers 15|] (noted in the text LKAG). 
Instead of calculating several magnetic configurations, this method, based on the magnetic 
force theorem (MFT)Q, H]) allows the evaluation of the MEI from one collinear configura- 
tion which is usually ferromagnetic. Computationally, this method is thus very attractive. 

Assuming infinitesimal rotation angles of the magnetic moments (limit of infinite magnon 
wavelength) is necessary to get the final LKAG formula for the MEI. However, one should 
note that this formalism is used for arbitrary big rotation angles (finite magnon wavelen gth' 



as well. Thus, many improvements of the formalism have be en p roposed recently: Bruno [18] 
proposed a renormalized MFT using the constrained DFT 19j leading to unrealistic high 
LDA Curie temperature (T c ) for fee Ni. The same effect has been observed using the proposal 
of Antropov[20]. Katsnelson and Lichtenstein proposed in their recent publication 21] a 

15] and the new renormalized theories 18|, |20] . They 



reconciliation between the old formalism 



have shown that the improvements proposed are well suited for the static response function 
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while the LKAG formalism is optimal for calculations of the magnon spectra. A more 



22 



-|28| 



rigorous approach is based on the calculation of dynamical transverse susceptibility 
which is computationally more involved. 

In the present contribution, we revisit the LKAG formalism and scrutinize one of the 
first assumptions assumed in the mapping procedure which has not been discussed yet. 
We demonstrate that an interesting issue occurs in the original mapping and thus in the 
majority of improvements as well. Avoiding the long wave or the infinitesimal rotation 
angle approximation, an anisotropy of the DFT MEI is obtained. This inconsistency is 
interpreted as a contribution to the DFT mapped part from high order MEI, such as the 
four-spin interactions, but behaving like bilinear terms. 

In our demonstration we follow the usual mapping procedure with three steps to consider: 
(i) definition of the classical Heisenberg model, (ii) evaluation of the DFT counterpart, (iii) 
mapping and extraction of the MEI. 

Classical Heisenberg model for pair interactions. As done in LKAG, we consider 
eq. [Hand determine the rotation energy of two spin moments at sites % and j, which are 
initially ferromagnetically aligned. Contrary to LKAG, here we assume different rotation 
angles for i and j. First, we determine the energy difference between this new magnetic 
state and the ferromagnetic one 



AE i+j = - J2 J in {e\ - 1) - J m3 {e) - 1) 

-JiM-ej-l) (2) 



where the z-axis refers to the quantization axis of the ferromagnetic environment and n 
and m to environmental atoms. Second, since we are interested in the MEI between atom 
i and atom j we subtract the interaction energies (AEi and AEj) of each atom with the 
environment. This is obtained after rotating only one of the two atoms, by the same angle 
as assumed for AE i+ j. 

AE . = _ Met - 1) - J fj (e? - 1) (3) 
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The final quantity which depends only on the MEI is thus given by 

AE (i>j) = AE l+j - AEi - AEj (4) 



AE (i,j) ~ ~ J ij 



1 + cos(0j) cos(9j) — cos(6 l j) — cos(9j) 



+ sin(6 l j) sia(9j) cos(fa — 



(5) 



if polar and azimuthal angles fa) and (6j, <pj) are introduced. In their work, LKAG cant 
the two spins by an equal angle 6 but in opposite directions i. e. by setting 9i = 9j = 9 when 
evaluating AEi and AEj while they cant the two spins by 9/2 and consider fa — <ftj = ir 



1 - cos(#) 



in agreement with 



when evaluating AEi + j. One then obtains AEuj) = —Jij 
LKAG. (Note that in the DFT counterpart expression LKAG use an angle 9/2 for AEi and 
AEj instead of 9). For small rotations 9 iy 9j eq. [5] simplifies to 

AE^j) m -Jij9i9j cos(fa - 4>j) (6) 

Magnetic pair interaction from DFT. This difference is directly given by 



dE(E - E F )An {hj) (E) 



dEAN(ij) (E), (7) 

with AN(i t j)(E) being the corresponding change of the integrated density of states (IDOS) 
and Ep being the Fermi energy. 

AN {hj) (E) = AN i+j (E) - AN t (E) - ANj(E), (8) 

Hence, ANi + j(E) is the change of the IDOS when both atoms i and j have their moments 
rotated. ANi(E) and ANj(E) are changes of the IDOS when only one moment is rotated. 
AN(ij}(E) is the change of the IDOS corresponding to the interaction energy between the 
moments i and j as expressed in eq. [5j 

Now, we can calculate every term in eq. [8] using multiple scattering theory and take 



advantage of the Lloyd's formula \2i 



AN(E) = --Im Tr nLs In (1 - G(E)AV), (9) 

7T 

where the trace Tr is taken over the site (n), orbital momentum (L) and spin (s) indices. 
Knowing the Green function G of the initial system describing the collinear magnetic state, 



this formula allows an exact determination of the change in the IDOS just by knowing the 
potential difference AV induced by the rotation of a magnetic moment. 

When rotating the magnetic moments of two atoms i and j, the interactive part of the 
integrated density of states according to eq. (jSj) is given by 



After taking the trace over n, the formulation giving the IDOS can be simplified into: 



which is equivalent to eq B.l from LKAG. Here we dropped out the argument E for reasons 
of clarity and the scattering t-matrices Ati and Atj describe all scattering processes at the 
isolated atoms i and j. At is defined by AV/(1 — GAV). 

The term GjiAtiGijAtj describes the scattering of an electron at a site j, the propagation 
to the site i from which it is scattered back to site j. It is a second order process which is 
expected to be very small compared to 1. A similar argument can be used for the denomina- 
tor. Indeed, if one makes a Taylor expansion of the denominator, terms like GAtGAtGAt 
would appear but are third order processes and thus are expected to be much smaller than 
1. 

After a first order expansion of eq. dO] we obtain 



The previous equation is expressed in the global spin frame of reference, i.e. the t- 
matrices have non-diagonal elements which is not the case of the magnetically collinear host 
Green function G. The MFT states that the spin-moment does not change upon rotation, 
meaning that the t-matrix within the local spin frame of reference of each atom does not 
change. Once calculated in the initial collinear state, the t-matrix is easily obtained: 





(10) 



ANfa) ~ -Im Tr Ls GjiAtiGijAtj 



(11) 



e bal (£) = \ [G{E)l + tlf f {E)U,. 



n 




(12) 



with U being a rotation matrix defined as following 



U 



cos(§)e-5* 
sin(f)e^ 



— sm(|)e 2 
cos(f)e^ 



(13) 
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and t l °£m and t di fj are equal to respectively t| ocal + tj° cal and t| ocal - fj° . From the 
t-matrix we subtract the initial one needed in eq. [TT] 



Atf^(E) = ^Af dlff (E) 



1 



cos(^) - 1 sin(^)e-^ 
sin(0;)e* - cos(^) + 1 



which is inserted in eq. [TT] leading to 



AN (hj) ~ -ImTrr 



47T 



(A + C)(cos(^)-l)(cos(^)-l) 



+2i? sin(6 l j) sin(^) cos(0j — 



after taking the trace over the spins with 



new 



(14) 



(15) 



A - GZj! At^ffG 3 ^ At diff , B - G^ At J diff GfAt diff , 



C - G l lAt diff G^At diff . 



(16) 



Mapping. Thus, the energy difference is given by 



AE, 



(ij) 



— J\{1 + cos(^) cos(#j) — cos(^j) — cos(9j)) 
— J 2 sin(6 l j) sin(^) cos(0. 



'3)1 



(17) 



where J\ = ^-Im TrL J Ef dE(A + C) and J 2 = TrL J F dE2B. This DFT expression 
is incompatible with expression flH]) calculated from the Heisenberg model, since two param- 
eters Ji and J 2 appear. Note that LKAG give only the expression for J 2 , which is also the 
expression used in the literature. However, it is only the correct expression for small angles 
9i, 6j, since J\ varies as 0f6j. We face here an important dilemma in determining the MEI, 
which, as we will show, results from higher spin interactions automatically included in the 
second order DFT approach. 

Let us evaluate the difference between the two terms: 



Ji - J 2 = y-Im Tr L 



dE 



(18) 



Since agreement with the Heisenberg model is only obtained, if J 2 = J\ or A + C = 2B, 
the difference J\ — J2 vanishes only if = G± i.e. for a non-magnetic reference system. 
This means that any magnetic system would lead to two possible values for the MEI. It is 




FIG. 1: The MEI J2 (circles) and J new = J2 + J2 ~^ Jl (triangles) calculated for bulk foe Ni (a) and 
bec Fe (c) with respect to the distance R with a being the lattice parameter. In the insets are 
plotted the terms A and B for Ni (b) and Fe (d). For reasons given below, the discrepancy between 
J2 and J ncw is much stronger for Fe compared to Ni. 

true that for magnetic excitations with tiny rotation angles or for what is called the long 
wavelength approximation (LWA), one gets rid off the first term in eq. [T7]but the error grows 
like (Ji — J2)(cos(#j) — l)(cos(9j) — 1). If the desired excited magnetic state is close to high 
values of the rotation angle then both terms J\ and J2 have to be considered. 

Using the full-potential Korringa-Kohn-Rostoker Green function method [32| within 
the local density approximation (LDA) 33] or the generalized gradient approximation 
(GGA)j34J, we evaluated these terms for usual bulk systems: Ni and Fe (see Fig. [[]) and 
found that J\ and J2 are, on one hand, relatively similar for Ni since it has very small 
magnetic moments (0.61 /is). On the other hand, Fe bulk is characterized by a stronger 
discrepancy due to its high bulk magnetic moments (2.3 hb)- 

In order to grasp some insight on the first term J\ we propose to consider from the model 
Hamiltonian side terms beyond the Heisenberg model which are expected to be implicitly 
included in the DFT counterpart. The additional terms can be obtained from a pertur- 



bation expansion of the Hubbard model jl, 



ail- 



The first terms which have been added 
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are the four-spin interactions, H 4 _ spin = - E m ^^? ^mr. w [(e m en) (epe ? ) + (e n e p )(e g e m ) + 
(e m e p )(e n e g )]/3. Calculating the energy difference (eq. [3]) due to the rotation of the atomic 
moments % and j leads to following further terms: 

-y (sin(^) sin(^) cos(& - </>,•)) (19) 

with K = J2 pH Kijpq- 

Obviously, one notices that the four-spin interactions with the uncanted environment 
spins behave for the ij pair like a bilinear term since only the moments i and j are canted. 
It is interesting to note that adding this term to the Heisenberg model brings an imbalance 
between the term proportional to the sin function and the one proportional to cos. We 
conclude that this mechanism is behind the observed anisotropy in eq. [T71 If we restrict 
ourself to the four-spin interactions only, the difference J\ — J2 would be given by 2K/3 that 
consequently would lead to a renormalization of from J2 to J2 + J2 ^ Jl that we represented 
as J new in Fig. [TJ This final result is without any doubt subject to modification as soon as 
higher order terms are included in the model Hamiltonian. The extraction of the exact MEI 
is thus a rather difficult task. As mentioned previously, since the moment of Fe is higher 
than the one of Ni, the discrepancy between the renormalized Jij and J2 is strongest for Fe 
(Fig. EE]) . 

We exemplify the effect of such corrections by evaluating the new Curie temperatures (T c ) 
by Monte-Carlo simulations. The extracted temperatures are not expected to be correct 
but are meant as illustrative examples for the effect of renormalizing the MEI. A major 
result shown in Table [I] is the large increase of T c with the new values of the MEI for 
Ni, Co and Fe. The difference between the old and new temperature gets stronger when 
increasing the magnetic moment of the host. Surprisingly, similar behaviors have been 



obtained by Katsnelson and Lichtenstein 



2 if ] when comparing the temperatures obtained 



usine; the renormalized method of Bruno [ljwith those of the old LKAG method. Obviously, 
the values obtained for Fe are too high and probably, one has to include higher order terms in 
the model Hamiltonian to lower T c . The values obtained with only J2 are probably sufficient 



for Fe due to a cancellation of errors that were described by Katsnelson and Lichtenstein 21]. 

By concluding we stress that the LKAG formula for describes correctly the MEI for 
small canting angles 9. In this case the spin-dependent t-matrices At of Eq. [T4l vary linearly 
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TABLE I: The Curie temperature (in K) for Ni, Fe and Co calculated with the LKAG formalism 
and with taking into account the 4-spin interactions. 



T c (K) Exp. .h j 2+ (^Ji) 

Ni(fcc-LDA) 631 374 458 

Co(fcc-GGA) 1388-1398 1520 1949 

Fe(bcc-LDA) 1045 1086 2062 

Fe(bcc-GGA) 1165 2791 



in 9, so that AE^j) is proportional to Jy 9 2 where 9 is an effective canting angle. All higher 
order interactions like K ijk i between 4 or 6 slightly canted spins therefore scale as 9 4 or 9 6 . 
As demonstrated, these Jy calculated by the LKAG formula include implicitly all multispin 
interactions of the canted moments with the uncanted environment atoms. It is for 
these reasons, that the calculated long-wave magnons and the spin stiffness constants agree 
very well with experiment. However, for larger transversal fluctuations of the moments 
the bilinear interaction J^- is no longer sufficient, and higher order spin interactions like 
the four spin interaction and the biquadratic coupling become important and have to be 
included explicitly in calculating T c and related thermodynamic properties. Since the spin 
splitting and At scale with the local moments M, these multispin interactions scale as M 4 
or higher and are thus more important for systems with large moments. In the paper, we 
have demonstrated the importantce of four spin interactions in T c -calculations for Fe, Co 
and Ni based on the LKAG formula for larger canting angles. 
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wishes to thank the Alexander von Humboldt Foundation for a Feodor Lynen Fellowship 
and D. L. Mills for discussions and hospitality. 
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